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Cubic invariants for two-dimensional degenerate Hamiltonian systems are considered by using variables of sepa- 
ration of the associated Stackel problems with quadratic integrals of motion. For the superintegrable Stackel systems 
the cubic invariant is shown to admit new algebro-geometric representation that is far more elementary than the all 
the known representations in physical variables. A complete list of all known systems on the plane which admit a 
cubic invariant is discussed. 



1 Introduction 

In 1935 Jules Drach applied direct method for search of the integrable Hamiltonian systems of two degrees 
of freedom, which admit a cubic second integral |2| . Recall, the direct approach leads to a complicated set 
of nonlinear equations, whose nonlinearity has no a priory restriction. We can attempt to solve these second 
invariant differential equations using various simplifying assumptions. The Drach ansatz for the Hamilton 
function H and for the second cubic invariant K 

H = p x py + U (x, y) , 

(1.1) 

K = 6w(x,y) (j^Py- Px ~^j ~ pl <Px,P yi x,y) 

yields ten new integrable systems. Recall, for the the fixed polynomial P(p x ,p y ,x,y) the potential U{x,y) 
and function w(x,y) be solution of some differential equations S. 
For the other natural Hamilton functions 

H=p 2 x +p 2 y + V(x,y) (1.2) 

the similar approach has been used by Fokas and Lagerstrom by Holt and by Thompson 0. Note, 
that the some Drach results have been rediscovered in these papers. 

The most complete classifications of known results was later brought together by Hietarinta in 1987 
Q . Recently |(| , two-dimensional hamiltonian systems with the cubic integrals were investigated using the 
Jacobi change of the time. In a complete list of all known systems was extended in comparison with 

The aim of this note is to study the Drach systems and some other degenerate systems on the plane with 
the cubic in momenta integrals of motion. We prove that eight Drach hamiltonians belong to the Stackel 
family of integrals Q and, moreover, seven of them are degenerate systems. 

Recall, the system is called superintegrable or degenerate if the Hamilton function H is in the involution 
with two integrals of motion I and K, such that 

{H,I} = {H,K} = 0, {I,K} = J(H,I,K). (1.3) 

Initial integrals and the constant of motion J(H, I, K) are generators of the polynomial associative algebra 
jl], 0, whose defining relations are polynomials of certain order in generators. 

Below we shall consider two-dimensional systems with two quadratic integrals of motion I\ = H, I2 = I 



and one qubic integral K. So, by the Bernard-Darboux theorem |17 the system with integrals I\, I2 belong 
to the Stackel family of integrable systems |@]. Therefore, let us begin with remaining of some necessary 
results about the Stackel systems H. 
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2 The Stackel systems 

The systems associated with the name of Stackel are holonomic systems on the phase space R 2n equipped 
with the canonical variables {pj,<?j}™ =1 - The nondegenerate n x n Stackel matrix S, with entries Skj de- 
pending only on qj 

detS^O, t^ = 0, j^m 
defines n functionally independent integrals of motion 

h = £ c ^ (p 2 j + u i) > c ^ = ' (2 ' 1} 

which are quadratic in momenta. Here C = [cjk] denotes inverse matrix to S and S J be cofactor of the 
element Skj ■ The common level surface of these integrals 

M a = {ze R 2 ™ : I k (z) = a k , k = l,...,n} 

is diffeomorphic to the n-dimensional real torus and one immediately gets 



*3 



dS N 2 



d qj 



= a i s n(li) ~ U j(h) ■ ( 2 - 2 ) 



i=l 



Here S(q\ . . . , q n ) is a reduced action function ||. For the rational entries of S and rational potentials Uj(qj) 
one gets 

s iiUto-e*) (23) 



where are constants of motion and functions (fj (qj ) depend on coordinate qj and numerical constants p"l| ] . 
The Riemann surfaces (2.3) are isomorphic to the canonical hyperelliptic curves 

k 

c j $ = II( A ~ e *) ' H = ( P(<lj)P3> ( 2 - 4 ) 

where the senior degree k of polynomial fixes the genus gj = [(k— l)/2] of the algebraic curve Cj. Considered 
together, these curves determine an n-dimensional Lagrangian submanifold in R 2 ™ 

C {n) : Ci(pi,qi) x C 2 {p2,q2) X ■■■ X C n (p n ,q n ) . 

The Abel transformation linearizes equations of motion on by using first kind abelian differentials on 
the corresponding algebraic curves JT^ ] . The basis of first kind abelian differentials is uniquely related to the 
Stackel matrix S 0,0. 

Now let us turn to the superintegrable or degenerate systems in the classical mechanics. One of the 
main examples of the two-dimensional superintegrable systems is the isotropic harmonic oscillator, which 
has many common properties with the Drach degenerate systems. Recall, for the oscillator the Hamilton 
function and the second integral of motion look like 

H =p\+pl+ql+ql, I =pl + ql-pl-q\. 
Obviously, the angular momentum 

K = q lP2 -p 1 q 2 = -( Pl — -—p 2 j . (2.5) 

is the third integral of motion. Two pairs of quadratic integrals I\ — H, I 2 — I and I\ = H, I 2 — K 2 are 
associated with the following Stackel matrices 

S=(\ \ ), and S=( r i a I), r 2 =x 2 +y 2 
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respectively. So, the corresponding equations of motion may be separated in the different curvilinear coor- 
dinate systems. 

For this degenerate Stackel system and for all other known superintegrable Stackel systems with quad- 
ratic integrals of motion the number degree of freedom n > g is always more then sum of genuses gj of the 
corresponding algebraic curves. In this case the number of independent first kind abelian differentials be 
insufficient for the inversion of the Abel-Jacobi map on C^ n \ 

To construct inversion of this map for the degenerate systems one has to complete a given basis of the 
differentials to the set of n differentials. We have some freedom in a choice of complimentary differentials 
and, therefore, we can associate the different Stackel matrices to one given Hamilton function |Q. By using 
first kind abelian differentials one gets superintegrable Stackel system with quadratic integrals only. Of 
course, we can try to add the second and third kind abelian differentials, but we do not know such examples. 

Below we prove that for all the known superintegrable systems with a qubic integral K the number 
degree of freedom n — 2 is more then sum g = gi + gi = of genuses gj — of the associated Riemann 



surfaces (2.4) too. The corresponding dynamics is splitting on two spheres 

Ci,2 : (J 2 = ai^X 2 + /3i, 2 A + 71,2 , #1,2 = 0, (2.6) 
where aj, (3j and jj be the constants of motion. 



In variables /zi. 2 (2.4) the additional cubic integral of motion for all the degenerate Drach systems looks 

like 

dot S / o> 2 dm \ 

K = Mi -77 jr • (2.7) 

S2is 2 2 \ at at J 



This generalized "angular momentum" gives rise to the first order integrals (2.5) or the third order poly- 
nomials in momenta depending on the Stackel matrices S and potentials Uj. In our case it will be qubic 
integral, which coincides with the initial Drach integral up to numerical factor. 

To consider nonlinear algebra of integrals of motion for the Drach systems we shall introduce new 
generators {N, a, at} instead of the two quadratic integrals 1\ = H , 1% = I, one qubic integrals K and the 



constant of motion J ( |l.3| ). Similar to oscillator these new generators have the following properties 

{N,a} = a, {AT, at} = -at, 

(2.8) 

{a, at} = $(/!, J 2 ), aat = *(h, J 2 ). 

Here generator iV(ii,/ 2 ), functions ^(Ji,^) and $(7i,/2) depend on the quadratic Stackel integrals only. 
Two other generators a and at be functions on the all three constants of motion I\, I2, K, such that 

K = p{h,h){a-a)), j(l l j 2i K)= < ^-. 



The relations ( |2.8| ) remind the deformed oscillator algebra, which is widely used for the superintegrable 
systems with quadratic integrals of motion However, instead of the usual quadratic algebra of integrals 

we shall get more complicated algebras of integrals. 

3 The Drach systems 

Let us reproduce corrected Drach results in his notations 

(a) U = — + f3x ri y r2 +jx r2 y ri , where r? + 3r, + 3 = , (3.1) 
xy j j 

P = (xp x - p y y f , w = x 2 y 2 /2 , 

U = + I 1 1SV±M (3 .2) 

^/xy (y - fix) 2 y/xy (y - ^x) 2 ' 

P = 3(xp x - p y y) 2 (p x + np y ) , w = xy(y - fix) , 



3 



(d) 



(e) 



(/) 



if)) 



(h) 



(k) 



U = 

P = 

u = 

p = 

u = 

p = 

u = 

p = 

u = 

p = 

u = 

p = 

u = 



a xy 



, 7 
(y — ax) 2 (y + ax) 2 



3(xp x - p y y)(pl - a 2 p 2 y ) , w = (y 2 - a 2 x 2 )/2 . 



\fyjx - a) y/yjx + a) Vx 2 - a 2 ' 

3p v [(xp x - PyV) 2 - a 2 pl] , w = -y(x 2 - a 2 ) , 



a P | 7 

fxy \fx Jv 



3p v Px {xpx - p v y) , w = -2xy , 

2x 2 + c jx 
axy + py = H — -^=^= , 

\J x 2 + c \J x 2 + c 

.2 „.„ ^ ... _ ^2 

n 



3p y (xPx - yp y ) , W = (x +c)/2. 



+ — 3mi) + 7(2/ — mx)(y — 9mx) , 



(y + 3ma;) 2 

(p x + impy) 2 (p x - 3mp y ) , w = -m(y + 3mx) , 



(» + -) 



mx,_ 2 / 3 



a + (3 (y — mx/3) + 7 (y 



I4mxy m 2 x 2 , 



2 WmpxPy m 2 p 2 



. rax 
w = -m(y+ — j 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(0 



P = 3p 2 x Py . 



U 



P 



(y-f)+/^- 1/2 + 7.*- 1/2 (y-H. 



3p x p 2 y + ppl, w = x. 



(3.10) 



Here a, (3, 7, /i, /?, a, c, and m be arbitrary parameters. In compare with Q we corrected function w in 
the case (g) (B/7) and revised potential U in the case (k) (3.9). Namely this corrected Hamiltonian is in the 
involution with the initial Drach cubic integral K (1.1). 

With an exception of three cases (a) ( |3.l|) , (h) (3^) and (1) ( |3.10| ), other Drach systems are degenerate 
or superintegrable Stackel systems. The separation variables associated with the pair of quadratic integrals 
{Ii = H, I2} are the Stackel variables. Equations of motion may be integrated in quadratures fn] , but these 
quadratures depend on the value of quadratic integral 12- Thus, instead of the solution of initial Drach 
problem related to integrals {H, K} we shall solve the associated problem with quadratic integrals {H, h}. 

In the case (h) (3.8) we also have the Stackel systems |u|. Only in this case (h) (3.8) dynamics is 
splitting on two tori and the number degrees of freedom is equal to the sum of genuses g = n = 2, such that 
the corresponding system is non-degenerate. 

In the case (1) ( 3.10D the Hamilton function coincides with the hamiltonian of the previous Stackel 
system (3.9) at p = 0. Here we shall not consider this generalized Stackel system at p 7^ 0. 

Below we shall consider the Drach integrals (1.1) up to linear transformations of the coordinates and a 
rescaling of these integrals. It allows us to remove some parameters in the Hamilton functions without loss 
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of generality. To associate the degenerate Drach hamiltonians with the Stackcl matrices S we can join these 
systems into the four pairs of the systems with a common matrices S. 



3.1 Case (a) 



In our previous paper MM, the first Drach system (3.1) has been related to the three-particle periodic Toda 
lattice in the center-of-mass frame. Namely, after canonical change of the time t = q n +i and the Hamiltonian 
H = Pn+i at the extended phase space 

dt = {xyy 1 -dt, H = xy ■ (H + 5) , 

and after further canonical transformation of other variables 

31 +'32 , . <?l+'92 11 -'92 , . , 91 -'92 

x = e 2 ; Px = ( Jpi - ip2 ) e 2 ; y = e 2 , Py = (pi + tp2)e 2 , 
the Hamilton function (|3.l[) becomes 

1 _ \/3 1 \/3 

ff = p? + p\ + {3e~2 qi ~ T q2 + ie ~2 qi + ~T q2 + Se^ + a . 

It is the Hamiltonian of the tree-particle periodic Toda lattice in the center-of-mass frame. The separation 
variables survive at the change of the time. Thus, for the first Drach system we can separate variables and 
integrate equations of motions in quadratures repeating the calculations for the Toda chain fllEfl . 
Later in || Thompson considered this system too. In fact, after point transformation 

e -i<t> e i4> 
^ = « 4 , p x = — — (p r ~ ip^r' 1 ) , y = re~ 10 , p y = — (p r + ip<j, r' 1 ) 



the Drach hamiltonian H (1.1) looks like 

2 

H = pl + ^ + U(r, ( j>), 



up to numerical factor. Namely this Hamilton function was studied in M and Q. The special substitution 
of the potential U(r, <j>) into the Drach equations leads to the following equation 

frM) = /w y" w , => f'f'- 2f"f -3/7 = 0, 

introduced in o. Of course, the same equation follows from the functional equation on the Toda potential 



3.2 Cases (b) and (e) 

Put \x = 1 in ([O). Let us introduce the Stackel matrix 



2 2 



1 1 



and take the following potentials 
(b) U x 



2a 



/3-2 7 



(e) 



U 1 = 2a + 2(l3 + 1 )q 1 , 



U 2 



U 2 



-2a 



/3 + 2 7 



-2a - 209 - 7)?2- 



(3.11) 



The corresponding Hamilton functions I\ (2.1) coincide with the Hamilton functions H for the Drach systems 
(3.2) and (3.5), after the following canonical point transformation 



{Qi ~ Q2) 2 



Pi -P2 

qi-q2' 



y 



(91 + 12? 



Py 



Pi +P2 

qi + q2 ' 
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The second integrals of motion J 2 (2.1) are second order polynomials in momenta. The third independent 
integrals of motion K are defined by (2.7), where 



{b) m = qipi, (i2 = 92P2, (e) fix = pi, (i% = p% . 

From the above definitions we can introduce generators of the nonlinear algebra of integrals ( |2.S| ) and verify 
properties of this algebra 



(b) 



N 



'2 



A-JH 



a=J + 4:VHK, (J = J-AVHK, 



an' = 16 {AH {(3 + 2 7 ) - (2a + I 2 ) 2 ) {AH (J3 - 2 7 ) - {2a - I 2 ) 2 ) 
{a, at} = -256VH (I 2 (I 2 - 2a) (Z 2 + 2a) - AH {f3I 2 - 4a 7 )) , 



and 



(e) 



2sf~H 



a = J + 2VHK, a) — J — 2\[H K , 



a a f = -16 (H {2a + h) - {/3 - 7 ) 2 ) (H {2a - h) + {(3 + 7 ) 2 ) , 
{a, a f } = -6AH 3 / 2 (7 2 H — (3 2 — 7 2 ) . 



3.3 Cases (c) and (g) 



Put a = 1 in ( p.3| ) and m = 1/3 in (3/7). Let us introduce the Stackel matrix 



S cff — 



1 1 

2 ~2 

1 1 



and take the following potentials 



G?) 



"9? 


7 


4 


2 ' 




a 

f — 


3 


91 



t r 2 



ag2 _ £_ 
A & 



U 2 



Aiql 



PQ2 



(3.12) 



(3.13) 



The corresponding Hamilton functions I\ (2.1) coincide with the Hamilton functions H (3.3) and (3.7), after 
the following canonical point transformation 



9i - 92 

X= , P x =Pl~P2, 



91+92 

V= 7, J Py=Pl+P2- 



The second integrals of motion I 2 (2.1) are the second order polynomials in momenta. The third independent 
integrals K are defined by (2.7), where 



(c) fa = qipi (i 2 = q 2 p 2 , {g) (11 = qipi , (1 2 = P2 ■ 
Generators and defining relations of the nonlinear algebra of integrals (2.8) look like 

h 



(c) 



N = 



a = J + 2^-aK, a [ =J-2y/-aK, 



a a f = {H 2 + AHI 2 + All ~ 4 "7) {H 2 - AH I 2 + All + ^ a P) , 

{a, a f } = -32VE {I 2 {H - 2h){H + 2J 2 ) - a/3(ff + 27 2 ) - a 7 (i? - 2J 2 )) 



G 



and 



(<?) 



N = 



h /3 



4 V 7 



a=J + 4jjK, a^ = J-4JjK, 



a T = 1/9 (8jH - I67J2 + 3/3 2 ) (3iJ 2 + 12iri 2 + 12lf + 16cry) 



{ a ,at} = 64Q)* 



7 \3/2 / (2J 2 + H)(4ryH -2AjI 2 + 3b 2 ) 1607 \ 
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3.4 Cases (d) and (f) 

Put a = 1 in (|3.4|) and c = 1 in (p]q). Let us introduce two the Stackel matrices 



( 1 1 ) 

1 1 

V a? 4 J 



9i 92 
1 1 

V if if y 



(3.14) 



and takes the following potentials 



(d) 



(/) 



2V2(a + £) 



9i 

7 (a + 2/?) 
2gi 4 



<7 2 = -2 7 - 



2V2(a - /3) 



'72 

7 (a -2/3) 
2 92 4 



The corresponding Hamilton functions I\ (2.1) coincide with the Hamilton functions H ( p.4| ) and (3.6) up 
to numerical factor, after the following explicit canonical transformations 



(d) 



q\ + ql (piqi ~ P2q2)qiq2 

X = — , p x = ■ 



2<M2 



? 2 

q{ - qi 



y = qi 92, Py 



,,s 91-92 2(pigi -p 2 92)^/91 92 . 

(/) a; = - = : — , y = y/qiq 2 , Py 



gigi + P292 
29192 

Pi 91 + P292 



2^/9192 



9i + 92 



V9192 



The second integrals of motion I 2 (2.1) are the quadratic polynomials in momenta. The third independent 
integrals K are defined by (2.7), where for the both systems one gets 



Hi = qipi 



H2 = 92P2 ■ 



Generators of the nonlinear algebra of integrals ( 2^8 ) are given by 

N=^/h, a = J + 2 v r hK, at = J — 2\fh, K , 
which have the following properties 

(d) a a) = 16 (h (27 -H) + 2(a + (3) 2 ) (l 2 (2 7 - H) - 2{a - [3) 2 ) , 

{a, a)} = {h (2 7 - H)(27 + #) + 2(a 2 + /3 2 )ff + 8a/3 7 ) , 



(/) 



1 



— (4/ 2 (a + 2/3) + (7 - 2iT) 2 ) (4/ 2 (a - 2/3) + (7 + 2#) 5 
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{a, a f } = -4V£T f (2/3 - a) (2/3 + a)I 2 + aH 2 + 2/3 7 iJ - 



cry 
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3.5 Cases (h) and (k) 

Put m = 3 in (|3.8|). Let us introduce two the Stackel matrices 





( qi 92 \ 




(<* 




Sh == 




, s.- 




T 




1-1 1 1 









and take the following potentials 



2 2 

(fc) £/ 1 =a + /3<? 1 + ^-, U 2 = -a + f3q 2 + ^. 



(3.15) 



The corresponding Hamilton functions I\ ( |2 . l| ) coincide with the Hamilton functions H ( [3.8] ) and (3.9) up 
to numerical factor, after the following explicit canonical transformations 

fi\ P2-P1 , (3gi+3g 2 ) 3 / 2 /3 Pl +p 2 

W + 54 + 16? P ^ \3 qi + 3 Q2 +V7fai-^ 



and 



P2 - Pi (3gi + 3q 2 ) 
4^7 54 



3/2 



i6 7 ' 



, P1+P2 

y/3qi + 3(72 



\/7(?l - 92), 



(fc) 



9i - 92 



Pi -P2, 



(9i + 92 ) 2 



Py = 



Pi + P2 

9i + 92 



Note, in the case (h) (|1|) we used non-point canonical transformation in contrast with other Drach systems. 

In the last case (k) ([}!]) integral of motion I 2 ( [2.l| ) is the second order polynomial in momenta. The 
third independent integral K is defined by (2.7), where 



/«i =pi and [i 2 =p 2 . 



Generators and defining relations of nonlinear algebra of integrals (2.8) look like 

h 



(k) N =^d=, a = J+^2^K, at = J-^p2F,K, 

a* = (2 7 (I 2 +a) + (H + P) 2 ) (2 7 (I 2 -a) + {H~ (3) 2 ) , 



{a, aT} = -47^^23 (H 2 + 2 7 / 2 + /3 2 ) 



In the case (h) ( |3.8| ) the second integral of motion J 2 (2.1) is the second order polynomial in momenta 
{pi,p 2 }. However, after the non-point transformation of variables this integral I 2 becomes the qubic in 
momenta { Px ,Py} Drach integral K (3.8). The corresponding dynamics is splitting on two tori and the 
third order polynomial (2.7) does not commute with the Hamilton function. Later this system has been 
rediscovered by Holt j^j . 



4 Other degenerate systems on the plane with a qubic integral of 
motion 

In this section we consider the Stackel systems on the plane with a qubic integral of motion defined by the 
following Hamilton function 

H = ^{pl+ P fj+V(x,y). 
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As above, the corresponding qubic integral will be written at the Drach form ([hi]). 

On the plane we know four orthogonal systems of coordinates: elliptic, parabolic, polar and cartesian. 
Thus, we reproduce all the known results [Q, [j) in correspondence with the type of the associated Stackel 
matrix QUI . 

The systems whose Hamilton functions separable in cartesian coordinates: 

{A) V = a(4x 2 + y 2 )+f3x + ^, (4.1) 

y Z 

P = PxP 2 y, W=|, (4-2) 

(B) V = a (x 2 +y 2 ) + ^ + ^, (4.3) 

u t \ x y 

P = {XPy-yp X )p x Py, W=—, 

(CO V = a(x 2 +y 2 )+(3 a XV , (4.4) 

(x z — y z ) z 

2 2 

x ~ ~ y 

P = (Xp y - y Px ) -P 2 y) , W= " , 

(D) V = a{9x 2 +y 2 ), (4.5) 

P = (x Py - Px y) p 2 y , w = ~Jg-> 
The systems whose Hamilton functions separable in parabolic coordinates: 

(F) V = (a+^— + —^)r- 1 , r = ^x 2 + y 2 , (4.6) 

r + x r — x 

P = (xp y -p x y)p x , w = _ jy> 

(G) V = (a+^f)r-\ (4.7) 

y 

P = (x py - p x y) 2 p x 

(H) V = (a + P VrTx + j V^x) r^ 1 , (4.8) 

P = (x p„ -p x y) i 2pi + 2 P 2 - JL_ - J. — ) , w = -- 



^2 

' v x , w = 

12 



(xp y ~p x y) (^2p 2 x + 2p 2 y - 



y/r + x \Jr — x j 6 
One system with the Hamilton function separable in polar coordinates: 

O V = P - 2 + - 1X ~ 5V , (4.9) 

y / x 2 +y 2 (dx + jy) 2 x 2 + y 2 {5 x + jy) 2 

(x -t- y ) (6 x -\- *y y\ 

P = (PxV - PyX)(jp x -Sp y ), W= — . 

It is new superintegrable system with a cubic integral of motion, which is a deformation of the degenerate 
kepler model. An application of the direct method ||, ||[ || or the Jacobi method || does not allows us to 
obtain this system (|4.9|). 
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Three exceptional systems whose qubic integral of motion K we can not rewrite in the "generalized 
angular momentum" form (|2.7|): 



(K) 

(L) 
(M) 



V 

V 
V 



w = ^Jxy . 



= a (y/x + 0y) , P = p 
= }'{4>)r- 2 , K = 



2(3 ' 

cos 4>p r — sin <p r _1 p^ + 



(4.10) 

(4.11) 
(4.12) 



(2/'(0) cos 4> - f{tj>) sin($) Pr + (3f (0) sin <f> + f{<j>) cos <j>) r~% . 



At the case (M) ( |4 . 1 2|) we used the standard polar coordinates {r,p r , <p,p$} and the function f{<f) has to 
satisfy the following equation 

/" (3/' sin <t> + / cos <t>) + 2f (2f cos <j) - / sin 0) = . 

At these exceptional cases the Hamilton functions (4.10 |4.11 ,4.12) are separable at the cartesian and polar 
coordinates, respectively. However, the Stackel potentials £7^2 are not polynomials in variables of separation. 

4.1 Cartesian coordinates, cases A-D 

Let us introduce the Stackel matrix 



Sa- 



A-D 



1 1 

2 2 

1 -1 



(4.13) 



and take the following potentials 

{A) 

(B) 

(C) 
(D) 



U 1 =8aq 2 1 +2(3q 1 , 



20 

Ui = 2aqf + 

Si 



U = 



aq\ 







U 2 = 2aqj + —j, 
12 

„ 27 
U 2 = 2aql + 

m, 

Mr 



(4.14) 



U 2 = 



aql 



Ui = 18a qf, 



U 2 = 2aq\. 



The corresponding Hamilton functions 1\ (2.1) coincide with the Hamilton functions H (4.1 4.3,4.5) and 
(0) if 

(A-B.D) x = q%, y = q 2 
or after the following canonical transformation 



(C) 



qi - 92 

X= , Px=Pl~P2, 



y 



qi + q 2 



Py = Pl +P2- 



The second integrals of motion I 2 (^J) are the second order polynomials in momenta. The third independent 
integrals K are calculated by (|2.7|), where variables 



(A) =pi [± 2 = q 2 p 2 , (B -C) /ii = gipi , M2 = P2 • 
determine the left hand side of the canonical algebraic curves ( |2.4| ). At the case (D) the variables 



1 n\ P 1 ^ 2 

{D) fii = p 2 qi — , \x 2 = p 2 q 2 



have not such natural algebro-geometric meaning. 
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Generators and defining relations of the nonlinear algebra of integrals (2J5) look like 

h 



{A- B) N 



such that 



(A) oo+ = 4(4a(2/ 2 + ff)+/3 2 )((2/ 2 -#) 2 -64a 7 ), 



{a, a" 1 } 
(B) a a f 
{a, a 1 } 

At the two last cases we have 
(C) N = 

a a 1 = 

{a, a*} = 



-128aV^ (2I 2 - H)(6I 2 + 64a 7 ) , 

2a 



f (2/ 2 + #) 2 - 64a/?) ((2/ 2 - H) 2 - 64a 7 
-16-/^2a (((2/ 2 - Hf - 64a 7 ) (2J 2 + if) 



+ ((2/ 2 + ff) 2 - 64a/?) (2/ 2 - i?) 
o= J + 2V-2a.fr, a f = J -2\/-2aK , 



2 V /= 2^ ' 

(2J 2 - ff) 2 - 2a/?) ((2/ 2 + H f + 2a/?) , 



-8V-2^((2/ 2 - H) 2 - 2a(3){2I 2 +H) 

+ ((2I 2 + H) 2 + 2a[3)(2I 2 -H)) , 



and 



(D) N 



6\/^2a" ' 

a a) = -4(2/2 - ff) 3 (2/ 2 + H) , 
{a,a f } = -96^/^2^ (2I 2 -H) 2 (4:I 2 +H). 



a^J + Gy/^^K, a) = J — 6V~2a K , 



At the case (D) (4.5) the quantum counterpart of this qubic deformed oscillator algebra has been used to 
study of the corresponding quantum superintegrable system [0 . 

4.2 Parabolic coordinates, cases F-H 

Let us introduce two Stackel matrices 

1 1 



S H = 



2 2 

91 92 



1 1 



(4.15) 



and take the following potentials 



a_ _p_ 
2 qi 2q* 



tj = _2 2_ 

J 2g 2 2q 2 ' 



(G) 
(H) 



u - -2L - JL 

1 2 91 4 ? r 



I7i = -4a-8V2/3gi, 



a (3 
U 2 = 4a + 8\/^2 7 g2- 



C/2 = _ ^ + ^l' 



(4.16) 
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The corresponding Hamilton functions I± ( |2.l| ) coincide with the Hamilton functions H (1.1, and (4.4) 
after the following canonical point transformations 



(F,G) 



X = 9l + q 2 , Px 



Pi 92 -_P2g2 
91 - 92 



2/ = 2\/-9l92 , 



(pi - p 2 )\J-q\q 2 



and 



(if) x = 9i + 92 j Px 



l glgi -p 2 g2 
2 9?"9 2 2 



, y = 2iq 1 q 2 , p y = - f 



9i - 92 



I Pi92 -P29i 



2 - qi 



The second integrals of motion I 2 (2.1) are the second order polynomials in momenta. The third independent 
integrals K are calculated by (2/7). where variables 



(F, G) fli = qipi (J, 2 = 92P2 , (H) fix = p\ , fx 2 = P2 



define the left hand side of the canonical algebraic curves (2.4). 

At all these cases (F,G) and (H) the generators and defining relations of the nonlinear algebra of integrals 



(2.8) look like 



(F-H) 



N 



2\ r H' 



a = J + 2VHK, at = J-2y/HK, 



such that 



^(8(3H-(2I 2 + a) 2 ) (87 H - (2I 2 - a) 2 ) , 



-2VH (l 2 {2h + a){21 2 -a)- 2H(p{2I 2 - a) + 7 (27 2 + a))) 



(F) aa f 
{a, at} 

(G) so* 

{a, at} = 2VH (I 2 (2/ 2 + a){21 2 - a) - 2a(3 H 



-—(4(lH+(2I 2 +a) 2 ) (4/3 H - (2I 2 - a) 2 



(H) aa f 



16rff(/ 2 -4a)+32 7 2 (if (7 2 +4a)-32/? 2 , 



{a, a 1 } = -64# 3/2 (#/ 2 -16/? 2 + 167 2 ). 



4.3 Polar coordinates, case I 

Let us introduce Stackel matrix 



1 



9r 2 1 



(4.17) 



and take the following potentials 



(I) 



U x =a + — 
9i 



U 2 = 



7 cos(g 2 ) - 5 sin(q 2 ) + p 
(5 cos(9 2 ) + 7 sin(g 2 ) 



(4.18) 



The corresponding Hamilton function I\ ( pTl| ) coincides with the Hamilton function H (4.E) if gi = r and 
q 2 = 4> are the standard polar coordinates on the plane. The second integrals of motion I 2 (2.1) are the 
second order polynomials in momenta. The third independent integrals K are calculated by (2.7), where 
variables 

(I) fj,x=px, P2 = P 2 (S cos(q 2 ) +7 sm(g 2 )J , 



define the canonical algebraic curves (2.4) 
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At S = 1 and 7 = the generators and defining relations of the nonlinear algebra of integrals (2.8) look 

like 



N 



a = J + 2J~T 2 K, a) = J-2J^T 2 K, 



a a ] 



[AH h - 4a/ 2 + )3 2 ) (4/ 2 2 -Apl 2 + 1), 
{a, at} = -8v / ^2((2/ 2 -p)/3 2 + (12/ 2 2 - 8p/ 2 + l)(i/-a) 



(4.19) 



This system does not contain in the list of the known integrable systems |4|, [| • At the cases (I) ( |4.9| ) and 
(M) (4.12) we have common leading part P of the qubic integrals K . However, at the case (M) we can not 



rewrite the qubic integral in the "generalised angular momentum" form. 



5 The Lax representation 

In (ll], |l2) we proposed some construction of the 2x2 Lax matrices for the Stackel systems with homogeneous 
Stackel matrices Jl^| and with uniform potentials Uj — U . The Drach-Stackel systems fall out from this 
subset of the Stackel systems. Nevertheless, we could construct the 2x2 Lax matrices for these systems by 
using various covering |ll| of the initial spheres C\p (p^). 

Here we consider the 4x4 Lax matrices for some Drach systems by using canonical transformations of 
the extended phase space, which induce transformations of the Lax matrices jl2|, Q. Recall, if the Stackel 
matrices S and S be distinguished the first row only, the corresponding Stackel systems are related by 
canonical change of the time q n +i = t and conjugated momenta p n +i = —H ]l2[| 



~ detS , 
at — - — — at . 



H 



~ H=^H, 



(5.1) 



det S det S 

Thus, starting with the Stackel systems related with matrix S = S cg ( |3.12 ) we can study systems associated 



with matrices S = Sf, e (3.11) and S = (3.15). Here subscripts mean the type of the Stackel matrices for 
the different Drach systems. 

The Stackel systems with the constant matrix S cg possess the following 4x4 Lax matrices |l6|, [llj 



£(A) 



Lx{\,px,q x ) 

L 2 (X,P2,q2) 

with independent 2x2 non-trivial blocks Lj{X). For instance, two standard blocks may be chosen 



(5.2) 



Pj 



A - qj 



MN 




Lj(X) 



V 



A 



A 



A 

mi 

A 



MN 



} 



Here <j>(\) is a parametric function on spectral parameter A and is the linear combinations of the Laurent 
projections 

According to [^2[ Q , canonical transformations of the extended phase space induce shift of the Lax 
matrices depending on the Hamilton function. Thus, by using one known Lax matrix C(X) ( |5.2] ) we can 
construct another Lax matrices. Namely, canonical transformations of the time (5.1) give rise the following 
shift of the corresponding Lax matrices 



£(A) = £(A) - H 



( 














\ 


a 





























V 








b 





/ 



a, b = ±1 or ± i , 



(5.3) 
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where values of the constants a, b depend on the chosen form of the blocks Lj{X). 

Below we present some Lax matrices constructed by designated above scheme. In the case (c) the Lax 
matrix is given by 

A 



£ C (A) 



Pi 



(A + ?i)(-f +79? A 2 



V 

so the spectral curve 
is a product 



-Pi 








ip2 



r(A, fx) 



i(A- 
det (C C (X) 



<?2 






«(A - <? 2 ) 
-ip2 



h 



" -y- /2 + — 



A^ 



= 
aA 2 



of the corresponding canonical Stackel curves (2.6). 
In the case (k) the Lax matrix is given by 

/ Pi A - qi 

7(A + <?i) 



V 



-a -pi 






W2 






i(A + g2 ) 
-ip2 



a 2 



7S =0. 



i/3 



/ 








°\ 


1 























V o 





i 





where H = I\ he the Hamilton function (3.9). As in the previous example the spectral curve 



is a product of the corresponding Stackel curves (£ 
In the case (b) the Lax matrix is given by 



7A 2 



h = , 



Pi 



pin 

A 

09 - 27)?r 



Pi<7i 



V 





P2g2 

A 

-p 2 + 08 + 2 7 )g- 





A 



P292 



f 








°\ 


1 























V 





1 





where H = 1\ be the Hamilton function (3.2). The corresponding spectral curve 



r(y, M ): det (a(A) - yl) 



= 



is a product 



h 



2a + I 2 (3 + 2 7 



A 



A 2 



2a 



0-27 



0. 



A A 2 

of the initial Stackel curves (2.2), which could be rewritten in the canonical form (2.6). 

All the s pect ral curves of these 4 x 4 L ax matrices £j , Ck and Cb give rise to the quadratic Stackel 
integrals I± t 2 (2T). The third integral K ( |2.7|) may be extracted from the same matrices by using multivari- 
able universal enveloping algebras fTl]] . In fact, this integral is a coefficient of the following multivariable 
polynomial 

trP 7r £(Ai) <g> £(A 2 ) ® £(A 3 ) ® £(A 4 ) , (5.4) 

where be permutation operator of auxiliary spaces corresponding to a Young diagram ir The 
nonlinear algebra of integrals (|]^) may be repr oduced by using the Poisson bracket relations between the 
Lax matrices ^C(Xj) [jTTfl . The formulae (5T) for the 256 x 256 matrices has been proved by using the 
computer algebra system Maple V. 
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6 Conclusion 



Let us discuss the list of all the known integrable natural Hamiltonian systems in the plane with a qubic 
integral |^|, || . We suppose that all these systems may be embedded into the family of the Stackel systems 
JLH , cither into the subset of the generalized Stackel systems [Q or these systems may be related to the 
Toda lattices and the Calogero-Moser systems [[l| As an example the last case (1) ( 3. 10[ ) of the Drach 
systems and the Fokas-Lagerstrom || model belong to the generalized Stackel systems 1 14|. The complete 
classification will be presented in the forthcoming publication. 

In this note we proved this proposition for all the Drach systems 0. Moreover, we rewrite the qubic 



integrals for the superintegrable Drach systems in common form (2.7). This generalized "angular momentum" 
may be used to construct another n-dimensional superintegrable Stackel systems with the cubic integrals of 
motion. For instance, let us consider the Hamilton function 



tt l t 2 , 2 , 2n , 7 + ^ . 1 



a (r 



(3(r + z) 



U(y/x) 



(6.5) 



where r = \/x 2 + y 2 + z 2 and U{y/x) be arbitrary function. The corresponding equations of motion are 
separable in the parabolic coordinates 



q 1 =r + z 1 q 2 = r- z, 
which related to the following Stackel matrix 



q 3 = arctan(y/a;) . 



S = 



110 
Ql 1 Q2 1 
9i~ 2 Q2 2 -4 



The Hamilton function (6.5) coincides with the Stackel integral I\ (2.1) if 

7 



Ui = — + - 



U 2 = 4 + - , U 3 = U(q 3 ) 

qi 92 



Thus we have integrable Stackel system with the independent integrals of motion 1\ 2 and ^3, which are 
quadratic polynomials in momenta. 

Canonical algebraic curves is defined in variables 



V\ = Pi qi ■ 



M2 = P2 92 



M3 = P3 



To substitute these variables in the "generalized angular momentum" (2.7) one gets additional qubic in 
momenta integral of motion K. In the initial physical variables this integral K looks like 



K 



(x +y )p z - 2zp z (p x x + p y y) - p z {p x x + p y y)" 



dH 



dH dH 



~0^{pxx + p y y)+Pz[p x +Py- x— - y— 



dx 



dy 



It will be interesting to understand the algebro-geometric origin of this "generalized angular momentum" 
(f 
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